Chaotic Rydberg atoms with broken time— reversal symmetry 
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The dynamics of Rydberg states of atomic hydrogen per- 
turbed simultaneously by a static electric field and a resonant 
microwave field of elliptical polarization is analysed in the 
quantum perturbative limit of small amplitudes. For some 
configurations, the secular motion (i.e. evolution in time of 
the elliptical electronic trajectory) may be classically predom- 
inantly chaotic. By changing the orientation of the static field 
with respect to the polarization of the microwave field, one 
can modify the global symmetries of the system and break 
any generalized time-reversal invariance. This has a dramatic 
effect on the statistical properties of the energy levels. 



One of the fundamental questions in quantum chaos 
is the correspondence between the classical dynamics of 
a given Hamiltonian dynamical system and the statis- 
tical properties of the energy spectrum of its quantum 
counterpart Universal fluctuation properties occur 
when the corresponding classical system is either inte- 
grable or chaotic. The spectral fluctuations for a classi- 
cally integrable system obey Poissonian statistics In 
another extreme, i.e. ergodicity, the statistical properties 
of quantum spectra can be modeled by ensembles of ran- 
dom matrices, namely with a GOE (Gaussian orthogonal 
ensemble) in the case of time-reversal invariant (or more 
generally, any other antiunitary invariance) systems or 
with GUE (Gaussian unitary ensemble) in the absence of 
such a symmetry For a generic Hamiltonian system, 
where typically chaotic and regular motions coexist, the 
fluctuations are expected to interpolate between these 
two limiting universal cases 1^^. While experimental 
observations of transition to chaos have been made for 
time-reversal invariant systems, the system studied be- 
low allows for the breaking of the antiunitary symmetry. 
As far as we know, this is the flrst experimentally real- 
izable example of such a situation in a quantum system 
(for wave chaos experiments see Q). 

Perturbed Coulomb systems are ideally suited for stud- 
ies of a manifestation of classical chaos in quantum me- 
chanics both theoretically and experimentally. A 
hydrogen-like atom in an uniform magnetic fleld reveals a 
smooth transition to chaos Q with increasing magnetic 
field. Similarly the motion of the electron in an atom 
driven by electromagnetic radiation undergoes a smooth 
transition to chaos with increasing field amplitude [Q. 
On the other hand, the anisotropic Kepler problem (mod- 



eling defects in semiconductors), becomes chaotic as soon 
as the effective mass becomes direction dependent [||. 
Different reactions to perturbations are just a manifes- 
tation of the inapplicability of the Kolmogorov-Arnold- 
Moser theorem ||] to a highly degenerate Coulomb prob- 
lem. 

In this letter, we consider the perturbation of the hy- 
drogen (H) atom by very weak static electric and reso- 
nant microwave fields. Classically, such fields may pro- 
duce chaotic dynamics in the secular motion of the atom 
(i.e. in the motion of the electronic ellipse). In the quan- 
tum mechanical language, states coming from a given 
manifold with fixed principal quantum number uq may 
significantly mix only with each other. Thus, we may ob- 
serve an "intramanifold" chaos, in which a finite number 
(possibly tIq) of quantum levels are involved (provided 
the coupling to the continuum is negligible). Intramani- 
fold chaos requires the strong mixing of at least 2 differ- 
ent degrees of freedom, here the total angular momentum 
and its z-component. 

The simplest system where such a behaviour has been 
expected is a hydrogen atom perturbed by uniform static 
electric and magnetic fields of arbitrary mutual orienta- 
tion [|jlO). At lowest non-vanishing order (small fields) 
the perturbation leaves a constant of the motion. By 
including second order terms, signatures of chaos have 
been observed in the quantum spectrum 10|. The nec- 
essary fields are, however, big and higher order terms - 
neglected in this approach - are of comparable impor- 
tance. Moreover, the high electric field leads then to fast 
field ionization which will blur the effects discussed. The 
situation discussed below does not suffer from this defi- 
ciency. Moreover it allows us to study the infiuence of an 
anti- unitary symmetry breaking on the level statistics. 

We consider a realistic three-dimensional H atom 
placed in a static electric field and driven by elliptically 
polarized microwaves. We define the z-axis as perpendic- 
ular to the plane of polarization of the microwave field. 
The Hamiltonian of the system in atomic units, neglect- 
ing relativistic effects, assuming infinite mass of the nu- 
cleus, and employing dipole approximation reads: 

H— ^ F{x cos tot + ay sin Lut) + 11 ■ r, (1) 

2 r 

where F, a and lu stand, respectively, for the ampli- 
tude, degree of elliptical polarization and frequency of 
microwave field while E denotes the static electric field 
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vector. The resonant driving of Rydberg states with 
principal quantum number no occurs when the atom is 
illuminated by a microwave field of frequency fulfilling 
ujQ — uj/uk = m, where m is an integer number and 
ujk — l/'^o frequency of the unperturbed Kepler 

motion. 

Because of time-periodicity, the Hamiltonian has a se- 
ries of quasi-energy levels, corresponding to the energy 
levels of the atomic system dressed by the microwave pho- 
tons. To calculate these quasi-energy levels, we employ 
a quantum perturbation method, whose details will be 
given elsewhere. In short, the method defines an effective 
Hamiltonian pl| in a given manifold, which takes into ac- 
count the direct coupling between the levels due to the 
presence of the static electric field and the indirect cou- 
pling through all levels of other manifolds, i.e. process of 
absorption and emission of microwave photons. The di- 
rect couplings gives a first order contribution of the static 
field while the indirect ones contribute at second order of 
the microwave field. The effective Hamiltonian inside 
a given hydrogenic manifold is thus the sum of a term 
proportional to E and a term proportional to F^. For 
arbitrary mutual orientations of the two fields, the two 
terms have incompatible symmetry properties and, when 
having comparable magnitudes, induce a global chaotic 
behaviour. The resulting matrix of dimension Hq is di- 
agonalized by standard routines. The method has been 
validated in the limiting case of combined weak static 
and linearly polarized microwave fields | p2[ . 

The semiclassical approach is an extension of the meth- 
ods used in our previous studies [ p^lp^ ] . It relies on quan- 
tization of the fast orbital motion of the electron keeping 
fixed the electronic ellipse. In the resonant case, this 
motion is pendulum-like. Taking the appropriate limit- 
ing case (valid for very weak fields) of the corresponding 
Mathieu equation solution yields the effective Hamil- 
tonian for the secular motion of the electronic ellipse, 
which is the classical equivalent of the effective quantum 
Hamiltonian discussed above. It expresses as: 



f2 

Tj _ _L0_p2 



Eq 



(2) 



where the constant — l/2ng term was omitted (i?scc de- 
notes the shift from the unperturbed energy level of the 
atom), and: 
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Fq = UqF, Eq = UqE, Lq = L/uq and Mq = M/uq are 
the scaled microwave and static electric field amplitudes, 
the scaled angular momentum and the projection of the 



angular momentum on the 2:-axis, respectively. The V' 
and (f> angles are canonically conjugate variables to the 
Lq and Mq momenta, respectively. Orientation of the 
static field vector with respect to the z-axis is given by 
the spherical angles p and 6. The expression for 7 looks 
complicated, but it is actually nothing but the compo- 
nent of the average atomic dipole on the static field axis. 
The explicit expression for Tmfi = r,„/ng is given in [p^ , 
Eq. (2.16). 

To calculate the quasi-energies of the system, one 
should perform quantization of the secular motion de- 
termined by the Hamiltonian (^). It has been done in 
simpler situations (e.g. H atom perturbed by linearly 
pol arized microwaves and a parallel static electric field) 
|l3| . Then the secular motion is integrable and its quan- 
tization straightforward. The present secular problem (2 
degrees of freedom) turns out to be non-integrable. 

Eq. (H) has some well defined scaling properties with 
the field strengths Fq and Eq. Let us define the reduced 
microwave strength T = FqUq/Eq = F'^tiq/E and the 
reduced Hamiltonian 



W — Flsccn^l Eq 
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The classical phase space structure depends only on the 
value of H and (beside the static field vector orienta- 
tion and the polarization of the microwave), but not on 
the detailed values of uq, E, F and the secular energy. 
Of course, weaker fields imply a slower secular motion, 
but this does not affect the structure of phase space. In 
the quantum mechanical picture, the energy splitting of 
a degenerate hydrogenic manifold also depends on abso- 
lute values of the fields, but the structure of the levels 
does not. 

For generic orientations of the fields, all symmetry 
properties are broken. It is only for ip equal to or tt/2 
that the system has some antiunitary invariance, under 
the combination of time-reversal transformation with re- 
fiection with respect to the YZ or XZ plane, see (|^). 
To observe the quantum signature of the breaking of the 
antiunitary symmetry, it is convenient to have a predom- 
inantely chaotic classical dynamics, as a transition from 
GOE to GUE statistics is expected. After rather exten- 
sive searches with various values of the parameters, we 
have found that the 2:1 resonance, i.e. ujq = m — 2, 
is more suitable than the principal one as it produces a 
more chaotic dynamics. The choice 9 w tt/A, a = 0.4 and 
J- — 5000 maximizes chaoticity of the system for reduced 
energy H in the range 8.5 — 9.5. A Poincare surface of 
section for these parameters (and p = 7r/4) is shown in 
Fig. |l| and displays chaos very predominantly. Moreover, 
the phase space structure does not change significantly if 
we set (p equal to or 7r/2 (for ip — the regular layer is 
slightly larger while for = 7r/2 it disappears). 

To obtain better statistics, the quantum levels 
were found for different orientation angles ip, 9 ^ 
{0.27r,0.257r,0.37r} and different no = 99 - 101 values. 
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keeping other parameters as in Fig. |l|. From each diago- 
naUzation, we took a fragment of spectrum correspond- 
ing to the quasi-energy interval 8.5 — 9.5 and unfolded 
the spectra to produce the statistical distributions. 

The cumulative NNS distribution is plotted in Fig. |^ 
together with the similar distributions corresponding to 
GOE and GUE §. It is apparent that the calculated 
numerical data traces between GOE and GUE, clearly 
demonstrating the breaking of the antiunitary symme- 
try. The data do not reach the GUE behavior because the 
corresponding classical dynamics is not entirely ergodic. 
To measure the departure from the GUE distribution, 
one may fit different theoretical distributions. A natu- 
ral choice is Berry-Robnik statistics |^ (for independent 
superposition of Poisson and GUE spectra) which allows 
one to estimate the relative measure q of the chaotic part 
of classical phase space. The best fit is for (7 = 0.94 which 
is in agreement with Fig. |l|. However, the agreement with 
the Berry-Robnik distribution is not perfect, see Fig. ||, 
which suggests that the deep semiclassical limit, required 
for Berry-Robnik statistics 1^, is not reached yet in our 
data. 

The statistical distribution proposed by Izrailev Q is 
more appropriate. Indeed, it fits the numerical NNS dis- 
tribution quite impressively (compare Fig. ||) yielding the 
"repulsion parameter" /3 = 1.47, roughly in the middle 
between the GOE (/? = 1) and GUE (/3 2) values. 

To focus on the long range correlations in the spectra, 
we study also the spectral rigidity, i.e. A3 statistics [Q. 
It gives an independent information about the relative 
measure q of the chaotic part of phase space. For a su- 
perposition of independent Poissonian and GUE spectra, 
one obtains |jl|,||: 

A3(L) = Af™'((l - q)L) + AG^^((?L). (5) 

Fitting Eq. (H) to the numerical spectral rigidity, results 
'mq — 0.94 which is in excellent agreement with the value 
of q obtained from NNS distribution. In Fig. ^, we also 
present the spectral rigidities obtained for antiunitary 
symmetry invariant cases, i.e. for <p = and = 7r/2. 
Their detailed analysis is left for a future publication. 
It is sufficient to say that, for = 0, the phase space 
structure is slightly more regular than one shown in Fig. |l| 
and consequently the fitted parameter q = 0.82 is smaller 
[for antiunitary invariant problem IS.^^^ is substituted 
by Af°^ in Eq. (|)]. On the other hand, for Lp = 7r/2, 
the classical phase space reveals more chaotic character 
and thus q = 0.95. Note that the saturation of A3, for 
all presented cases, appears for L w 40 Hj. 

We have here considered small but finite field ampli- 
tudes to stay well within the applicability range of the 
quantum perturbation theory. We expect, however, a 
similar behavior for larger amplitudes. While we present 
numerical data for ng « 100 for better statistics, a clear 
signature is also observed for no ~ 50. There, the ex- 
pected mean level spacing is of the order of few MHz (for 
Fg R:! 5 X 10~^) making the experimental resolution of 
the spectrum feasible. 
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FIG. 1. Poincare surface of section (at <j) = 0) of the clas- 
sical secular motion, Eq. for Ti. = 8.5, with parameters 
of the fields T = 5000, uo = 2 (i.e. m = 2), a = 0.4. Ori- 
entation of the static electric field vector ip — 9 — 7r/4. Note 
that, for the parameters chosen, not the whole {Lo,tp) space 
is accessible. 

FIG. 2. Cumulative level spacing distributions W{s). In 
panel (a), the solid line represents the numerical data for 
no = 99 - 101, = 5000, c^o = 2, a = 0.4 and different static 
electric field vector orientations ip, 9 = {0.27r, 0.257r, O.Svr} 
(there are about 6000 spacings); dashed line (hardly visible 
behind the solid line): the best fitting Izrailev distribution 
with its parameter (5 = 1.47; dotted and dash-dotted lines: 
GOE and GUE distributions respectively. Panel (b) shows 
a fine-scale representation of the deviation of the numerical 
level spacing distribution from the best fitting Izrailev distri- 
bution in terms of the U{W{s)) — U {Wizraiiev{s)) vs. W{s); 
the transformation U{W) = arccos \/l — W is used in order 
to have uniform statistical error over the plot [5]. The up- 
per and lower noisy curves represent one standard deviation 
from the calculated numerical data which thus lie in the mid- 
dle of the band. The solid curve is the best Berry-Robnik 
distribution, q = 0.94. 

FIG. 3. Spectral rigidity, A3 compared with random en- 
sembles predictions. Panel (a) shows the results for numeri- 
cal data (solid line) in the broken anti-unitary invariance case 
(with the same parameters as in Fig. ^) together with the fit of 
Eq. ^ (dots, barely visible behind the solid line, fitted value 
q = 0.94). Dashed lines indicate Poisson, GOE and GUE pre- 
dictions as indicated in the figure. Panel (b) shows numerical 
data for an antiunitary invariant case with the same param- 
eters as in Fig. ^ but for Lp — and no = 97 — 102 for the 
upper case and = 7r/2 for the lower line. The fitted values 
for the fraction of chaotic phase space volume are q = 0.82 
and q — 0.95, respectively. 
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